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Abstract. The reaction matrix of a cavity with attached waveguides connects 
scattering properties to properties of a corresponding closed billiard for which the 
waveguides are cut off by straight walls. On the one hand this matrix is directly 
related to the S-matrix, on the other hand it can be expressed by a spectral sum 
over all cigcnfunctions of the closed system. However, in the physically relevant 
situation where these eigenfunctions vanish on the impenetrable boundaries of the 
closed billiard, the spectral sum for the reaction matrix, as it was used before, fails 
to converge and does not reliably reproduce the scattering properties. We derive here 
a convergent representation of the reaction matrix in terms of eigenmodes satisfying 
Dirichlet boundary conditions and demonstrate its validity in the rectangular and the 
Sinai billiards. 

Recently, there has been some interest in the application of the reaction-matrix theory 
of Wigner and Eisenbud []T| and the projection-operator formalism of Feshbach 0, 
originally developed for the description of nuclear collisions, to chaotic cavities with 
attached scattering channels || £|, ||, ||, [7j . Such models are frequently used as paradigms 
of chaotic scattering || [|, |10|, [TTJ and found important experimental realizations by 



electron transport through open quantum dots u2, 13|, lasing optical micro-cavities 



14] , [L5| , |16f and scattering of microwaves in resonators with attached waveguides |L7|, ||. 
It is known that the quantum scattering in the open system shows signatures of the 
classical dynamics in the closed system. For example, conductance fluctuations of open 
quantum dots are different for systems whose closed counterparts have integrable, fully 



chaotic or mixed phase space [0, [13| . In fact many statistical results for quantum 
chaotic scattering rely on this connection as for their derivation an ad-hoc formulation 
of the scattering problem in terms of an effective non-Hermitian Hamiltonian is used 



2Tf , which describes the dynamics inside a closed system and an additional coupling 



of the corresponding eigenstates to some scattering channels. Classical chaos enters via 
the random-matrix assumption for the Hermitian part of this Hamiltonian [ITR |HS|, 



However, while the effective Hamiltonian appears naturally within the formalism from 
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nuclear physics, it is not a priori clear when this formalism applies to some given billiard 
system and how the parameters of the two different models are related. 

Therefore it is important and interesting to develop a thorough understanding of the 
connection between the scattering properties of billiards, in particular the S-matrix, and 
the properties of the corresponding closed system, i. e. spectrum and eigenfunctions. For 
example, today it is well known that the spectrum can be found from a secular equation 
involving the S-matrix ^|, [10| . For scattering from the outside of convex billiards in 
1Z 2 this scattering approach to quantization allows even for a mathematically rigorous 
formulation |22| . 

Unfortunately, the opposite direction is less profoundly understood. Here, the 
unitary S-matrix is related to a Hermitian reaction matrix, and this can in turn be 
expressed as a sum over the internal spectrum with coefficients reflecting the behavior 
of the internal wavefunction at the boundary separating billiard and waveguide. From 
the physical point of view, and in particular for the aforementioned applications 
fl2| , |T3| , |H] , 15], [16 , 17| , the most natural situation is a wavefunction which vanishes 



outside the billiard and on the boundary (Dirichlet b. a): Electrons in a quantum dot 
are depleted from the boundary by a high negative gate voltage, the radiation field is 
restricted to the optical cavity by total internal reflection or additional mirrors, and in 
microwave resonators the metallic walls do not admit the electrical field. However, for 
this of all choices of boundary conditions no consistent representation of the reaction 
matrix as a spectral sum is known. We derived a formal expression |3|, |[] but found 
both numerically and from semiclassical estimates that it fails to converge. Recently, 
this conclusion was confirmed with different methods |7]]. 

It is certainly possible to circumvent this problem: One option is to ignore the 
divergence and to restrict the spectral sum by hand to some set of levels. Naturally this 
must fail in a generic situation, but it can give reasonable results in some special cases, 
e. g. when the spectrum has a doublet or band structure [0]. A second option relies on 
the fact that in principle all boundary conditions providing a self-adjoint Hamiltonian 
are admissible for defining the closed system, see j7| for a nice explanation of this point. 
In particular Neumann b. c. with finite wave function but zero derivative at the interface 
were considered previously 0, J|, || || [7j. In this case the spectral sum converges, and the 
resulting reaction matrix successfully reproduces numerical or experimental scattering 
data. Nevertheless it would be very unsatisfying, could the formulation of a proper 
reaction-matrix theory for billiards not be based on the physically relevant boundary 
conditions. 

To remedy this situation we establish in the present paper a convergent expansion 
of the reaction matrix for Dirichlet billiards with attached waveguides. Our main result 
is Eq. ( |To| ) below, where the coupling constants entering this series are given. Before 
we get to this equation, we recall some results from previous work. Following Eq. ([1"5|) 
we give a derivation of this formula, discuss some interesting aspects related to it and 
show with two examples how it works. 

We consider the same "frying-pan" setup as in Refs. || £§, i. e. a cavity with an 
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Figure 1. A scattering system consisting of a cavity and an attached waveguide of 
width b is shown with bold lines. The coordinate system is chosen such that x is the 
parallel and y the transversal coordinate with respect to the waveguide. 



attached waveguide as shown in Fig. [I]. The wave function must vanish on the boundary 
of the system. To obtain the corresponding closed system, Dirichlet boundary conditions 
are imposed also on the line x = 0. It is assumed that this point of separation is well 
inside the waveguide such that for — a < x < oo the wavef unction can be expanded in 
transversal modes 



<t> x (y) = y/2/bmn(Xiry/b). (1) 

A condition on the minimum value of a will be given later. 

Energy eigenf unctions in a billiard satisfy the Helmholtz equation (A + k 2 )ip(r) = 



with k = \j2mE/h 2 and r = {x,y}. In particular, the scattering states ip x are in the 
channel region given in terms of the S-matrix as 

V*(r, k) = £ ( Sx , t Ae -i*v- + Sy, x e +ik ^) , (2) 



X' 



where k\ = yk 2 — (Xir/b) 2 is the longitudinal momentum. There is a finite number 
A = [kb/ir] of open modes with real k x . All other modes are evanescent and give 
contributions to the scattering states which decay exponentially into the waveguide. 
It is known that their influence is negligible unless the energy nearly coincides with a 



threshold for the opening of a new channel [10, o]. For simplicity we will disregard 
evanescent modes here. Then the A x A S-matrix is unitary, and therefore it is related 
to a Hermitian reaction matrix |] K by 

K= -i(I + S)-\l-S) (3) 
S = (I + \K)-\I -iK) . (4) 

As billiards are invariant under time-reversal, S and K are symmetric matrices. 

On the other hand we consider the eigenfunctions of the closed billiard which vanish 
at x = and can therefore be represented as 

*n(r) = ^2u n> x<px{y)sia(k n ,xx)/k n> x (-a < x < 0) . (5) 

A 

The corresponding eigenvalue is k n , and k n ^\ follows in analogy to k x . Again because of 
time-reversal symmetry the coefficients u U)X can be chosen real. These coefficients are 

§ Despite some minor differences in the definition we use this term here because of the strong analogy 
to |ffl|. 
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obtained by projecting the normal derivative of the eigenfunction at the interface on 
the transversal modes 

u n>x = f dy<t> x {y) [d/dx^ n {r)] x=0 . (6) 

J 

It is our goal to represent at arbitrary wave number k the reaction matrix K (and 
thus via Eq. (|3|) also the S-matrix) in terms of the discrete set {k n ,u ntX }. For this 
purpose we consider the Green's function for the closed billiard which is defined by 
(A + k 2 ) G(r; r', k) = S(r — r'), as this quantity can at the same time be represented in 
terms of the S-matrix and by a spectral sum. Indeed we have [Q 

G M k)= Zh]p*M (7) 

n=l K ft n 

= s A(r, k) gxx{k) ^ y {r', k) {-a < x' < x < 0) . (8) 

A, A' 

with 



and 



s A (r, k) = <p x (y) sin(k x x) / J k x (9) 



g(k) = (I + S(k))- 1 . (10) 

The S-matrix of the scattering system is implicitly contained in Eq. (0), both in the 
scattering states ip\ and in the matrix g(k). We can extract this information from the 
normal derivative of Green's function at the interface, projected to transversal modes. 
After inserting Eqs. (fj) and (pi) into Eq. (Rl) we have 







2 



lim \im ——G XjX/ (x,x') = iJk x kx>y29\\" (-$\",\' + S\",\>) (H) 

x'-*-0x-+-Qd x O x > V ^Jf 

Using Eqs. fliTf ) and (^|) we see that the r.h.s. reduces essentially to an element of the 
reaction matrix 

1 d 2 
K xx , = —== lim lim — — G x ,\'(x,x). (12) 

Because of the symmetry of Green's function with respect to its spatial arguments the 
result does not depend on the order of limits, just the singularity occurring for x = x' 
must be avoided. Comparing Eq. (|l^) to the spectral decomposition of Green's function 
Eq. (0) one naively expects that the reaction matrix has a spectral representation in the 
form 

K ^ = ^22 h 2 h 2 ( 13 ) 

n=l K K n 

with coupling constants 



W n . x = u n)X /^itk x . (14) 

However, this is not the case, as it was found in Refs. 0, [7| that with Eq. ( |T4"D the sum 
(Pf) fails to converge. Convergence is lost because the individual terms in the (absolutely 
convergent) series Eq. (|7|) grow too much upon the differentiation required by Eq. ([T2|) . 
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We will show in the following that Eq. fll^D is still valid, albeit with modified coupling 
constants 

Wn,A = l/- i V ,A- (15) 

Keeping in mind that differentiation and summation do not commute in Eq. (|7]), we need 
a different strategy which allows to make use of this series anyway: We will represent 
the second derivative of Green's function in Eq. ( jl"2|) by the value of Green's function at 
some shifted point in space. To implement this idea we return to Eq. ([TT|) and bring it 
by some formal manipulations into the shape of Eq. ([8]), projected onto the transversal 
modes. Then we get a representation of the reaction matrix which is equivalent to 
Eq. ( |T2D but does not require differentiation, namely 

Kw = i gx\" (—8\",x' + Sx",y) 



X" 



kxkx> £ ^ 1 9xx»^= (Vve+-/ 2 + S»',ye-^) 



'k x k X ' Gx,x> (~n/2k x , -n/2k x ,) . (16) 

Now we can safely use the spectral decomposition of Green's function Eq. (|7D together 
with Eq. (|) to obtain 

1 . fk n: xTT\u ni x . f k n> X' 7r\ U Ul X' 



' = ? W^i 8m [—2 ) Sm (172 j k^, ' (17) 

This is indeed equivalent to Eq. (|T3|) with the coupling constants given in Eq. (]T5|) . 

Some comments are in order at this stage. The point X\ = —n/2k\ in Eq. (|16|) 
corresponds to the first maximum of the partial wave A. Although no differentiation was 
involved, the normal derivatives of the wave functions u n> x appear in Eq. ( |TTD because 
they determine the amplitudes of the partial waves at these maxima according to Eq. fl5|) . 
However, the full reaction matrix K can be obtained by the outlined procedure only, if 
the points x\ are inside the waveguide for all A = 1, . . . A. Otherwise Eq. ([TBI) breaks 
down. This implies a restriction which is at its strongest for the minimum of k\ at 
A = A. We can represent the wave number as k = (A + jb with < k < 1 and 
find k\ = \/2Ak + k 2 t\ jb. Hence, our expression for the reaction matrix will be valid 
provided that 

b/2a < \/2Ak + k 2 . (18) 

This condition will always be violated at the threshold energies for the opening of new 
scattering channels where k = 0. On the other hand, if we avoid these singular points by 
fixing k to some positive value, the condition will always be satisfied in the semiclassical 
limit A — >• oo. In some sense these restrictions on the validity of our approach are similar 
to those allowing to neglect evanescent modes || [Tl], |JJ. From Eq. ( |T8"D we conclude that 
the situation which is most favourable for reaction-matrix theory in Dirichlet billiards 
is a closed billiard which extends very far into the waveguide b <C a. Intuitively this 
should be clear, because then the shape of the billiard resembles a scattering system. 
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A second comment concerns the uniqueness of the suggested procedure. According 
to Eq. and because of the special geometry we consider, any point —a < x < oo 
inside the waveguide region can be used to extract information about the S-matrix, while 
we have selected only those points where some partial wave A has maximum amplitude. 
This may appear like an arbitrary choice. Indeed, if we sample Green's function at 
sufficiently many points, arbitrarily placed inside the waveguide, it is in principle possible 
to compute the S-matrix or the reaction matrix. However, the resulting expression will 
not have the canonical form of Eq. (|l3| ) and no coupling constants for individual levels 
can be defined in this case. Therefore the connection to the formalism of Refs. [|], |2|, [yj 
would be lost within such an approach and consequently it would not be very useful. 
In other words, the main accomplishment of our theory is not the mere possibility to 
compute the S-matrix from data obtained in the closed billiard, rather it is the fact that 
the resulting expression is still of the form of Eq. (|i~3|). 

Further we note that the coupling constants in Eq. (|T5|) reduce to the naively 
anticipated Eq. fll~4]) for states n which are close to the energy shell of scattering, i. e. 



W n j ~ u ni i/\/k\ir for k n ~ k. This explains why Eq. ([[4]) works (up to some degree 
of accuracy) for a cluster of almost degenerate states 0, and why it also provides a 
consistent answer for the mean coupling strength in the semiclassical regime [[|: The 
approximation k n ~ k was always made from the outset. In general this approximation 
is not justified and Eq. ( |T4"D breaks down. 

The crucial difference between Eqs. fll~4] ) and ([15]) is their behavior for n — > oo which 
decides the issue of convergence. To discuss this question we need the semiclassical 
estimate for the average behavior of the boundary functions (|w n ,A| 2 ) An,A ~ 4/A which 
we derived in Eqs. (34), (38) of Ref. || for chaotic billiards (see [^3| for more general and 
accurate estimates of this type). A denotes here the area of the billiard. Substitution 
into Eq. (|l|) leads to 

— — — k n ~ n I (n->oo), (19) 

i. e. Eq. (^) is indeed absolutely convergent. 

We find it illuminating to apply presented theory to a simple toy model, where 
all relevant quantities are known in closed form. Consider a rectangular billiard with 
side lengths a and b which yields a half-infinite empty waveguide of width b if one of 
the walls is removed. This is an integrable system since the transversal modes of the 
waveguide are decoupled. Therefore we can restrict attention to some particular mode 
A. The corresponding diagonal element of the S-matrix is found to be 

S x = -e 2ikxa . (20) 

This is to be reproduced by the spectral sum Eq. (0). The normalized eigenf unctions 
of the closed billiard are 

y) = 2/ Vabsm(fi7cx/a) sm(Xny/b) (21) 
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Figure 2. In the energy interval with 10 open modes the numerically exact S-matrix 
is compared to the result of Eq. (|l3| ) for the quarter Sinai billiard shown in the inset 
(b = 7T, R = 7r/2, R + a = 57r/4). (a) Shown are the eigenphases 9\(k) (A = 1, . . . 10) as 
a function of the wave number k. Full line: exact result. Dots: approximation based on 
all 1,000 eigenfunctions up to fc max = 36.27. Vertical line: Range of validity of Eq. ( |l3| ) 
as given by Eq. (|l^). (b) Shown is the deviation of the total phase 9 = J2\ from 
the exact result for some fixed wave number k — 10.50 (vertical line) but with varying 
cutoff fc max for the included billiard eigenfunctions. 



and therefore we have the longitudinal momentum k^x = fin /a, the boundary function 
u^x = d'Z/akft and the coupling constants 

». \—™(tf\- ( 22 ) 




2akx / 

Substitution into Eq. (|I3|) yields after some straightforward transformations 

akx ~ l-cos(^f ) 

k x = —yJ2 ( t 1 l ak \ - (23) 
7T 2 {akx/n) 2 - /i 2 

According to the two terms in the numerator we split this expression into two series 
which can separately be evaluated with the help of Eq. (1.445.6) from Ref. pi] . Upon 
recombination of the two results we have 

K x = cot k x a , (24) 

which is precisely the reaction matrix corresponding to the S-matrix Eq. (p0|) via Eq. (Q). 

Now we test the outcome of Eq. ([13]) in the Sinai billiard. S-matrix and eigenstates 
were computed with the methods described in Ref. [0. For Fig. |2| we have chosen 
a large number of open modes, A = 10, and a geometry which does not allow for 
narrow resonances. This is a very demanding regime, as a high number of states is 
expected to contribute. Nevertheless, in Fig. ^a we show that a reasonable agreement 
can already be obtained from the 1,000 lowest internal modes. The agreement improves 
as we increase the number of states (Fig. |b). For fc max = 150, corresponding to 15,575 
modes, hardly any deviation is visible on the scale of Fig. [|a (not shown). In contrast, 



Reaction matrix for Dirichlet billiards 



8 



from Eq. (|HD we get garbage irrespective of fc max (also not shown). Note that the error 
in an individual eigenphase 9\ of the S-matrix is at its largest when 6\ = 0. These 
are the eigenvalues for a billiard with identical geometry but Neumann b. c. at x = 0, 
while at 9\ = ±7r the Dirichlet billiard is quantized and we get perfect agreement. Had 
we chosen Neumann b. c. as starting point for the evaluation of the reaction matrix, 
the result would show relatively large deviations close to 6\ = ±7T. In this sense the 
different boundary conditions are complementary. Finally we remark that from Eq. (fL8|) 
we cannot expect agreement in the interval 10 < k < 10.025 in Fig. §a. Indeed we 
observe in this region a few points which are very far off the exact result. 

Summarizing we have shown how and under what circumstances the reaction matrix 
of a cavity with attached waveguides can be represented as a spectral sum of a closed 
billiard with the physically relevant Dirichlet boundary conditions. 
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